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In discrete-coordinate quantum models the kinematical observable of position need 

not necessarily be chosen local (i.e., diagonal). Its smearing is selected in the nearest- 

Ph| neighbor form of a real asymmetric (i.e., cryptohermitian) tridiagonal matrix Q. Via 

Gauss-Hermite illustrative example we show how such an option restricts the class 

a 



of admissible dynamical observables (sampled here just by the Hamiltonian) . 

1 Introduction 



> 

Virtually any textbook on quantum mechanics pays attention to the point particle 
moving along straight line. Tacitly [TJ, it is assumed that the argument x of its 
wave function ip(x) G L 2 (R) coincides with an eigenvalue of an operator of the 
^ . (by assumption, observable) spatial position. The wave functions ip(x) carry the 

standard probability-density interpretation and describe the system in the so called 
x— representation. 

For Hamiltonians H = p 2 + V(x) in which the external potential is local, such 
an approach is natural. Due to the second-order differential-operator form of H, one 
easily determines the energies. The calculations may be further facilitated by the 
transition to approximations. Typically [2], the equidistant N— point Runge-Kutta 
discretization 

xeR — ► Xj = xf K) = xf K) + h ■ j , h>0, j = 0,l,...,iV-l 

converts operator H into a real and symmetric tridiagonal N by N matrix. 

For non-local external potentials the concept of the coordinate loses its guiding 
role of a link between the quantum and classical pictures while the kinetic energy 
operator p 2 may still be simplified via Fourier transformation J 7 : ip(x) — ¥ ip{p)- 
Thus, unless one has to deal with some other observables defined as functions of x, 
the computing costs of the mapping T remain acceptable. 

The feasibility of the calculations worsens for the more complicated Hamiltonians 
H 7^ p 2 + V so that the role of the quality of the approximations increases. In 
particular, the Runge-Kutta grid points Xj = x, may prove far from optimal in 
practice [3]. Many non- Runge-Kutta discretizations are being proposed and used in 
physics and quantum chemistry, therefore [3]. 



Table 1: The sample of the closed- form grid points x\ at N < 5 
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In our forthcoming methodical considerations one of the most typical samples of 
the non-Runge-Kutta one-dimensional grid points will be selected in the form of the 
iV— plet of the zeros of the N— th Hermite polynomial [5], 

x eR — ► Xj = xf P) , j = 0.1, ...,N-1. (1) 

In a way explained in Appendix A we felt puzzled by the conflict between the amend- 
ment of the numerical efficiency gained by the transition from x, to x, (cf. 

Refs. [5j |6] for details) and the loss of the closed formulae for Xj (with a few 
low— N exceptions sampled in Tabled]). A softening of the latter disadvantage (i.e., 
of the purely numerical character of the amended values x, at higher N) will be 
proposed here, therefore. 

In section El first of all, the numerical values of x, will be reinterpreted as the 
implicitly defined eigenvalues of a suitable tridiagonal matrix Q with elementary 
matrix elements. It will be argued there that one of versions of such a next-to- 
diagonal form of Q^ is in fact often used during the evaluation of values Xj in 
numerical practice. 

In section |3] (complemented also by Appendix B) our specific choice of manifestly 
non-Hermitian N by N matrices Q^ will be advocated. Their extreme simplicity 
will be identified there with the main criterion of applicability of the increasingly 
popular use of representations of quantum observables (sampled here by the position 
and energy) in the so called cryptohermitian picture. 

In section HI the main formal aspects and consequences of the use of the grid x, 
in such an implicit representation will be explained, in detail, via the first nontrivial 
N = 4 model. The abstract picture of kinematics (represented by the matrix Q ) 
and dynamics (represented by the related admissible families of Hamiltonians H^ N ' 
and/or of the so called metric operators Q^ N >) will be given there the concrete forms 
in which the construction of matrices Q^ in the simplest (viz., tridiagonal and 
pentadiagonal) sparse-matrix forms will be paid particular attention. 

In the final section |5] the resulting change of perspective transferring emphasis 
from the simplicity of quantum kinematics to a potentially better balance between 
the simplicity of kinematics and dynamics will be summarized. 



2 Formalism: Kinematics 



(HP) 



2.1 An implicit definition of grid points x- 

Let us consider the classical orthogonal Hermite polynomials 

H {x) = l, H x (x) = 2x, H 2 {x) = 4x 2 -2, H 3 {x) = 8 x 3 - 12 x 



(2) 



and make our present grid points x- unique by defining them as the roots of 
equation 

' 1 ^ P) )=0, J = 0,1,. 



h n I ^ x j 



iV-1. 



(3) 



Naturally, such a definition of the grid points is purely numerical, at the larger 
lattice-sizes N at least. At the finite N, there is in fact no practical necessity of 
insisting on the multiplicative-operator nature (i.e., on the diagonal matrix form and 
representation) of the position-operator Q = Q^ N \ Thus, once we re-read our implicit 
lattice-definition (J3J) as a vanishing-condition for the determinant of recurrences for 
Hermite polynomials we may immediately replace the traditional and, in essence, 
purely numerical diagonal representation of the operator Q^ N ' by its following, well 
known tridiagonal-matrix alternative 
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(4) 



The proof of the coincidence of the eigenvalues of this matrix with the roots of Eq. ([3]) 
is easy - it suffices to recall the recurrences for Hermite polynomials j7]. Subsequently, 
we may even specify the related eigenvectors of matrix (J5|), i.e., the smeared position 
eigenstates in closed form, 
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(5) 



In such an overall setting one has to circumvent several conceptual as well as purely 
technical obstacles. Some of them will be discussed in what follows. 



2.2 A trivial Hermitization of Q^ 
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As long as we have Q^ ^ Q^ we must abandon the "friendly" real vector space 
1{( F ) = M. N as unphysical. At the same time we may introduce a diagonal N by N 



matrix tt with positive matrix elements along its diagonal, 
This enables us to define matrix 
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(6) 



possessing the same eigenvalues as Q. As long as the original position operator Q 
itself is a tridiagonal matrix, we are allowed to require that our grid-points-preserving 
similarity transformation ((6]) has the Hermitization property, 
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(7) 



Formally speaking [8] we may now treat the new matrix q as acting in an idealized, 
"physical" Hilbert space "H which remains isomorphic to WL N . 



3 Formalism: Dynamics 

3.1 Hamiltonians 

The standard textbook picture of quantum dynamics may be perceived as living in 

( p) 
Hilbert space 7i . In this space the time evolution may be assumed generated by an 



arbitrary "effective" Hermitian Hamiltonian matrix f)o 
backward map 
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will then define the non-Hermitian Hamiltonian acting in the old "friendly" space 
TiS F '. The original Hermiticity of f) 0e is strictly equivalent to the Dieudonne's j9] 
quasi-Hermiticity relation rewritten in the double-dagger-superscripted notation and 
in terms of the abbreviation Vt\flQ := @ , 
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b;; x ! . (9) 

Under the implicit methodical assumption that the two operators of observables q 
(= position) and \y Q J (= Hamiltonian) in V, are "prohibitively complicated", the 
determination of at least some of the properties of the quantum system in question 
may still prove simpler in the unphysical, auxiliary Hilbert space T-L^ where the 
isospectral operators of observables Q^ (= position) and HqJ (= Hamiltonian) 
appear manifestly non-Hermitian. Thus, their direct use in computations (e.g., of 
their spectra) may happen to remain well motivated [TO] . 



3.2 The introduction of the third Hilbert space 

Let us now assume that at any fixed integer a > 1 and in all of the text of paragraphs 
12.21 and 13.11 one replaces the diagonal (i.e., trivial) matrix Of, by a fully general 
matrix Q a (to be called a Dyson's map jlO]) such that the product Q a = £l a £l a 
(i.e., a certain less trivial metric) remains sparse and strictly (2a + 1)— diagonal. 
Under this assumption the kinematics of our schematic quantum system may still 
be represented by the same position matrix Q( N ) defined in %^ F \ Naturally, the 
same specification of kinematics (i.e., of the grid points) also reappears, in the new 
physical Hilbert space H a , as carried by the new, manifestly Hermitian matrix q„ . 

Naturally, the transition from a = to a > will allow us to change the dy- 
namics in general. This is the key idea of our present paper. Thus, the choice of 
any a— dependent dynamical input information represented, say, by some "effective" 
Hamiltonian \y a J remains fully at our disposal. In our present exemplification of 
the theory this Hamiltonian will be chosen as a matrix which is Hermitian in Ti^ 
and which determines the observable properties of our hypothetical, a— parametrized 
quantum system. 

The related generalized, a— parametrized version of the non-diagonal Dyson's 
map Q, a will be again, in general, non-unitary, implying that Q a ^ I. This means 
that in a close parallel to the preceding two paragraphs where we used a = 0, we 
shall be forced to declare the "friendly" Hilbert space Hq "false" and unphysical. 

In such a case, in a way outlined in [8], it proves useful to introduce the third, 

(s) 
"standardized" Hilbert space Ha . By definition, the latter space may coincide with 

its "friendly" partner H^ as a vector space of kets \ip). At the same time it will 

certainly differ from it by its Q — dependent definition of the Hermitian conjugation 

of operators. 

Formally, we may merely replace the zero subscripts in Eqs. flBJ), flS]) and flU]) by 

their generalized forms, o — >a, while calling the positive definite operators Q a = 

Q a the "Hilbert-space metric operators" [10J. Indeed, by construction, the physics 

described in the third Hilbert space H a remains strictly the same as in H a . The 

unitary equivalence between Hilbert spaces n> a and W a may be also emphasized 

by the notation recommended in Ref. [8] and defining the linear functionals (a. k. 

a. the Dirac's bra- vectors (ip]^ or rather, in the compactified notation of Ref. [5], 

"brabravectors" ) in the less usual Hilbert space H a via the Hermitian-conjugation 

operation which is a— dependent and reads 

T Q (S) :|^> -> (^ s ):=<(^| = (^|e a . (10) 

The values of the inner products ((ip\4>) in the "standardized" physical Hilbert space 
H^ may be also evaluated inside the auxiliary Hilbert space W^ F ' since 

«V#) !*.*(*> = (V>i©l0)L^) • (ii) 

Thus, we never leave the standard quantum mechanics. The l-ft 8 '— specifying metric 
must, of course, have all of the required properties (i.e., basically, 6 = 6^ > in 
the present scenario with N < oo [10] )■ In other words, the requirement of the 
consistency of the theory may be expressed, inside the auxiliary Hilbert space 7i^ F \ 
as the compatibility of our matrices of observables (viz., position Q = Q( N > and 
Hamiltonian H = H a J) with the Dieudonne's constraints 

&® = Q® (12) 



H^e = HQ 



(13) 



They represent just the necessary conditions of the (hidden) Hermiticity of the re- 
spective operators of observables. 



4 An illustrative example with N = 4 

4.1 General non-diagonal isocoordinate mappings Q^ 

At one of the simplest nontrivial choices of N = 4 the present "smeared" (i.e., non- 
diagonal) operator of position reads 
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The related exhaustive (i.e., four-parametric) solution of Eq. f JT2|) is then easily 
found to read 







(4) 
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d 1/2 p 1/2 d+ 1/8// l/12p+l/48fc 



(14) 



Once we select fc = 1 and // = p = <i = 0we re-obtain the above-discussed a = 
metric O with elements [0o]n = 1, [©0)22 = 1/2, [@o]33 = 1/8 and [60)44 = 1/48. 
In an opposite extreme of the fully general four-parametric metric ( TT4l) , one can find 
the complete set of the underlying isocoordinate mappings Vl^ via the factorization 
of 9 = fttft. 

The inspection of the second independent constraint (TTBl reveals that in the 
trivial diagonal-metric example with O = O an important subset of the admissible 
Hamiltonians H will be composed of arbitrary diagonal real matrices. In such a 
setting one might speak about the (discrete) wave functions in energy-representation. 

4.2 Tridiagonal one-parametric metric 

At k = 1 and p — d = the a = 1 metric found from Eqs. (1T41) or (TT2l) reads 
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(15) 



It might be factorized, in particular, into products O = Q^Q using the lower-two- 
diagonal ansatz for the real factor Q and evaluating its matrix elements in the re- 
current manner starting from the lower corner. 



Once we wish to avoid the recurrent factorizations (which do not lead to any 
nice formulae at general N), a direct algebraic method may be used whenever the 
off-diagonal part of the metric remains small. For illustration let us notice that in 
such an approach one reproduces the metric of Eq. (TT5|) by the approximate N = 4 
and a = 1 matrix 



n(n) 



which is exact up to the fourth-order uncertainty factor 1 + 0(fi A ). In the same 
perturbative spirit we may accept the larger uncertainty 1+0 (/i 2 ) in the approximate 
inverse 
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and obtain the first-order correction to formula (J7J) above, 
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(16) 



In this precision the Hermitized position matrix is tridiagonal. The grid-point eigen- 
values evaluated from the approximate position matrix ( ITBT) coincide with the exact 
values at the small /x (cf. Figured]). 
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Figure 1: The quality of approximation of the /x— independent spectrum of Q^ A ' by 
the /x— dependent grid-point eigenvalues of the Hermitian position matrix q^ 4 \fi) 
considered in the leading-order tridiagonal small— /x approximation ffTB"]) . 

The consistency of the perturbation calculations with respect to the small /x may 
be reconfirmed by the inspection of Figure |2] in which one sees that the smallest 
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eigenvalue of the exact one-parametric a = 1 (i.e., tridiagonal) metric f JT5|) ceases to 
be positive at the not too large values of |/i| ~ 0.3. The exact formula 



eS 4) (/i)>0 iff /IE (-IM>,(M>) 



/io = 1/^6 + 2 V6 « 0.3029054464 
is also available due to the exact solvability of the model at iV = 4. 
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Figure 2: The \i— dependence of the quadruplet of the eigenvalues of the tridiagonal 
candidate ffT5l) for the metric. 



4.3 Two-parametric pentadiagonal metric 

Our above commentary accompanying Figure [2] means that at the larger /i the re- 
quirement of the positivity of the metric can only be satisfied at a > 2, i.e., in a less 
trivial physical domain T> of the available variable parameters k, /i, p and d entering 
the fully general formula of Eq. ( fT5|) for the metric. 
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Figure 3: The boundaries of the domain of positivity of the pentadiagonal candidate 
for the metric. 

The simplicity of our N = 4 model enables us to make this argument more 
quantitative. The explicit a = 2 choice of the two variable parameters fi and p in 



the metric 
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can be used for this purpose. In a preliminary graphical analysis we reveal (cf. Figure 
[3]) that the two-dimensional version of the physical domain T> has the piecewise linear 
boundaries. A more rigorous explanation can be also provided since the boundary 
lines may be identified with the loci of zeros of the secular polynomial S(/x,p, 77) = 
det(6;> (/i,p) — tj) at r\ — 0. Once we evaluate 

S(/i, p, 0) = det 9^ (/i,p) = 

= tL ~ i/SpV - 1/6 p 3 + l/16p 2 + 1/12 p 4 + l/48p- i-/i 2 + i-/i 4 
768 o4 64 

we reveal that the polynomial S(fi,p, 0) factorizes. This explains the existence of the 
two left-right symmetric pairs of the straight nodal lines which cross at fi — 0. The 
upper one is prescribed by the equation 



H'upper r ^~ l -£\P P upper) -"-IP 



2^-4 



(in Fig. [3j this is the upper pair crossing at the positive, maximal admissible p uvver 
1.112372435). Similarly, the lower pair of lines given by the equation 



Pic 



±(P - Plower) = ± P + 



2^ + 4 



crosses at the minimal, negative pi OW er = —0.1123724357. 

Along the line p = 0, Fig. [3] correctly reproduces the position of the two bound- 
aries (fI7|) of the interval of the positivity of all of the eigenvalues of G^ (/i) as 
displayed in Fig. [2J The same picture also explains why this interval of positivity of 
the metric grows with p > up to p m aximai = 1- 



5 Summary 

In accord with all textbooks on quantum mechanics the set of the operators of 
observables (sampled here by the position q and energy f)) must be kept Hermitian 
in a certain physical Hilbert space of states (here, in H^). In such a setting the 
starting idea of our present paper was that such a Hilbert space may be replaced by 
any unitarily equivalent alternative (in our notation, by H,^). 

The majority of the existing applications of such an idea is built upon the unitary, 
Fourier-like correspondence J 7 between bra or ket elements of the two spaces. In our 
present paper we rather made use of an indirect, manifestly non-unitary correspon- 
dence Q involving just the ket vectors of the two spaces [8j. As far as we know (cf. 
also several additional comments collected in Appendix B), the first use of such a 
correspondence (called Dyson map) found its applications in nuclear physics [TU] . 



Our present form of Hilbert spaces %^ or 7i^ has been chosen finite-dimensional, 
generated via a systematic iV— point discretization of the real line of coordinates. The 
underlying idea was that in such a setting the efficiency of the calculations need not 
necessarily be optimal when we diagonalize, by brute force, all of the matrices in 
question (i.e., here, the pair of matrices q and f)). We imagined that even the Her- 
miticity of the matrices of observables need not be a condition sine qua non. We 
argued that the non-diagonal forms of these matrices, whenever sufficiently com- 
pact, would offer a way towards a minimization of the time needed for the numerical 
evaluations of phenomenological predictions. 

For illustrative purposes we started from the choice of the "kinematic" grid-point 
eigenvalues x, of a matrix q. We imagined that the latter matrix may be further 
simplified via a nonunitary Fourier-like transformation f2 -1 , 

n (P) _> n W . (i 8 ) 

In the latter, naively chosen Hilbert space the isospectral analogue Q of the position 

appears non-Hermitian but the remedy is easy. In the subsequent "Hermitization" 

step 

n W __> n (S) (19) 

one simply endows the old space with a new, nontrivial metric 8^7, i.e., with an 
alternative, amended inner product. 

Our present attention has been paid to the definition of the third space 1-L^ S \ 
We showed that and how the trick (TT8l) plus (TT9l) opens a new paradigm in the con- 
structive analysis of quantum dynamics by making the pullback Q — y q internally 
consistent. We showed that even though our pre-selected operators q or Q of the 
coordinate were not diagonal, their respective Hermiticity inside w- p > and w- s ' has 
still been, constructively, guaranteed. 

On such a purely kinematical background the alternative dynamics were assumed 
given by the variable class of energy observables f) . We explained how these Hamil- 
tonian matrices may be related to their cryptohermitian avatars H, and how their 
eigenfunctions might be assigned the standard probabilistic interpretation. Thus, in 
contrast to the existing literature as reviewed, say, in [TTJ and starting from the wave 
functions in x— representation, we did not postulate any Hamiltonian in advance. We 
merely fixed its eligible class compatible with the pre-selected, smeared-coordinate 
representation of the wave functions. 

Our main emphasis has been put on a minimal smearing of the position. For 
the entirely pragmatic purposes we selected the N— plet of grid points Xj adapted to 
the optimal precision of the (discrete, approximate) integrations. Our matrices Q^ 
have been chosen tridiagonal, representing just one of the most elementary lattice 
grids given as the Hermite-polynomial zeros. We expect that such a strategy will 
optimize the precision of approximants, especially when evaluating or prescribing the 
norms and/or inner products of wave functions. 

Summarizing, we proposed here the kinematical choice of the cryptohermitian 
position Q( N ^ accompanied by the dynamics-determining specification of one of the 
eligible metrics Q a . In such a framework we indicated the possibility of varying the 
dynamics and/or of moving towards a systematic and exhaustive construction and 
classification of all of the alternative, a— numbered physical Hilbert spaces, be it 
Ha and/or its unitary equivalent, less usual version Hi with nontrivial metrics 
given, in our illustrative example, in closed form. 
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Appendix A: A comment on choices of grid points 

One of the roots of the popularity of the representation of the quantum-state ket- 
vectors \ip) in their wave-function (alias "coordinate-representation") form ip(x) G 
L 2 (R,fi) may be seen in the user-friendliness of the formula giving their norm, 



il)*(x)ip(x)uj(x)dx . (20) 

In a way similar to the definition of the Hilbert-space inner product 

/•oo 

(ip\4>) := / ip*(x)(f)(x)u(x)dx 



these integrals contain a weighting function u(x) which is most often set equal to a 
constant or a power of x or a Gaussian. 

The necessity of the evaluation of these integrals is ubiquitous. Its precision 
represents one of the limiting factors of the practical applicability of the theory. 
In particular, various empirical tests are very often passed by the so called Gauss- 
Hermite quadrature [5] 

poo N-l 

/ e~ x2 f(x) dx =Y, Wj f{xj) + error[N, f^ 2N \0] • (21) 

In such a Gaussian-weighted approximative formula the size of the error term is very 
often found acceptably small over the relevant class of the integrands. Its grid-point 
quantities Xj are identified with the the roots of the well known Hermite polynomials 
(this identification is sampled by our present Eq. ([3]) above). The necessary weight 
factors are also obtainable in compact form, Wj — Wj = C(N)/[H(N — l,Xj)] 2 . 
The more detailed information (i.e., the proportionality factor C(N), etc) may be 
found not only in standard monographs (12] but also, directly and free of charge, via 
web H31II3. 

Appendix B: A comment on the concept of crypto- 
hermiticity 

In the Dieudonne's paper [9] operators Q satisfying relation f JT2|) at a suitable positive 
and self-adjoint "metric" were studied and given the name of "quasi-Hermitian" 
operators, Q G C^ H '. In their full generality and without additional assumptions, 
they were emphasized to exhibit peculiar and strongly counterintuitive, truly "non- 
Hermitian" properties. In particular, once they were not assumed to be Riesz op- 
erators or compact symmetrizable operators, their spectrum was shown not to be 
necessarily real. Moreover, their adjoints were shown non- quasi-Hermitian in gen- 
eral, Q ] i C {QH \ 
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In the historical perspective it may be considered rather unfortunate that more 
than thirty years later, Scholtz, Geyer and Hahne [TO] decided to choose, incidentally, 
the same name for the well-behaved and very narrow compact-operator subset C^ SGH ^ 
of the Dieudonne's quasi- Hermit ian set. This is the reason why we use here the 
terminology inspired by Smilga [15] and why we call the elements of the set C^ SGH ^ 
cryptohermitian operators. 

The most natural and physics-motivated extension of the cryptohermitian class 
£{sgh) jg ^e Mostafazadeh's [16] set of the so called "pseudo-Hermitian" operators 
£( PH ) f or which, roughly speaking, all of the invertible and self-adjoint operators 
6 compatible with Eq. f lT2|) (and denoted, therefore, by another symbol, say, V) 
may remain indefinite. A potentially equivalent name of VT— symmetric operators 
is usually reserved, by its proponents [TTJ, [TOJ US], to the narrower class of models 
where V coincides with the operator of parity and where it may be interpreted as 
the indefinite metric in Krein space [20] . 

One of the most compact approaches to the applications of VT— symmetric op- 
erators in quantum theory has been reviewed by Carl Bender |llj . He considers 
merely the intersection of the cryptohermitian and VT— symmetric classes of oper- 
ators of potential observables while picking up a unique metric in the product form 
q(cb) _ -pQ r-j-^g k e y mer jt f this approach is that the factor C with the property 
C 2 = I may be interpreted as a charge of the system. The alternative name for 
C as proposed in [21] (viz., "quasiparity" ) has not been given the meaning of an 
observable and found much less users in the literature, therefore [22J. 

Acknowledgement 

Work supported by the GACR grant Nr. P203/11/1433, by the MSMT "Doppler In- 
stitute" project Nr. LC06002 and by the Institutional Research Plan AV0Z10480505. 



13 



